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SUMMARY 

Thrust losses due to the inability of dissociated combustion gs.ses 
to recombine in exhaust nozzles are of primary interest for evaluating 
the performance of hypersonic ramjets. Some results for the expansion 
of hydrogen-air combustion products are described. Combustion air was 
preheated up to 3300° R to simulate high-Mach-number flight conditions. 
Static- temperature measurements using the line reversal method and wall 
static pressures were used to indicate the state of the gas during 
expansion. 

Results indicated substantial departure from the shifting equilib- 
rium curve beginning slightly downstream of the nozzle throat at stagna- 
tion pressures of 1.7 and 3.6 atmospheres. The results are compared with 
an approximate method for determining a freezing point using an overall 
rate equation for the oxidation of hydrogen. 


INTRODUCTION 

The high combustion temperatures encountered in jet engines for 
high-Mach-number flight lead to extensive dissociation of the combustion 
products. Since residence time in the exhaust nozzle is short, the com- 
bustion products - because of necessarily finite reaction rates - cannot 
follow the shifting equilibrium during expansion. If these recombination 
rates are slow, large thrust losses can occur. These losses are shown 
for a stoichiometric hydrogen-air system in figure 1 for the limiting 
cases of frozen and equilibrium expansions (ref. l) . The losses can be 
appreciably greater for hydrocarbons or systems involving condensable 
exhaust products . 

A number of mathematical solutions have been recently applied to 
the expansion through a nozzle with reaction (refs. 2 to 6) . These are 
essentially stepwise integral methods and thus far have involved only a 
single reaction. For fast computers, no limitations are seen to exist 



for extending the methods to simultaneous reaction rates. The difficulty- 
in applying these solutions is the present lack of knowledge of reaction 
rates, and relaxation rates in general, at the conditions of interest. 
Reaction- rate data obtained by measurement of concentration gradients in 
flat flames (refs. 7 to 12) and in shock tubes (e.g., ref. 13) may be 
applicable to the nozzle problem. The solutions must be checked, how- 
ever, against actual measurements in a nozzle before any confidence can 
be established. The assumption that rates determined at conditions far 
removed from equilibrium can be applied to conditions near equilibrium 
leaves some doubts, since perturbations from a Maxwellian distribution 
may occur during radical recombination reactions (ref. 14). 

Only a few experiments have been conducted to study recombination 
in a nozzle. In one of these experiments, where inlet conditions could 
be varied to produce flows varying from near-equilibrium to near-frozen 
flow, measurements were made for the recombir ation 2N0 2 ** N2°4 ( re * > * 6 ) • 
These showed good agreement with a complete solution of the flow. 

The present work describes measurements of static pressures and 
temperatures for the hydrogen-air system in a convergent -divergent ex- 
haust nozzle and is an extension of previous work described in reference 
15. The operating conditions are designed tc approximate ramjet flight 
at a Mach number of 6 at an altitude of about 100,000 feet. Fuel-air 
ratios up to 1.4 times stoichiometric were si udied. The nozzle is coni- 
cal with an exit-to-throat area ratio of 9.5:1 and a throat diameter of 
3.08 inches. 

The calibration procedure used for the j ine reversal pyrometer is 
described by Donald R. Buchele in appendix B, entitled "Calibration of 
Pyrometer with Arc Source. " 


M 


TEST FACILITY 

Air from the laboratory system was meteied through a standard ASME 
orifice, heated, and passed through the test section, where fuel was in- 
troduced, burned, and the combustion products expanded through a super- 
sonic nozzle. The flow was then cooled and c ucted to the laboratory 
exhaust system. The facility is shown schematically in figure 2. 

The storage- type heat exchanger (ref. lb) used to preheat the air 
consisted of a packed bed, 4 feet in diameter and 10.4 feet deep, insu- 
lated with refractory and insulating brick. The bed was packed with 
alumina balls 3/8 inch in diameter. It was 1 eated by burning gasoline 
and air and passing the combustion products c.own through the bed. The 
bed was banked continuously at 1200° to 1500 c F when not used for test- 
ing. During a test, the bed was heated to the desired operating temper- 
ature, after which the airflow through the bed was reversed. The heater 
will supply air at 3.6 atmospheres at a temperature up to 3300° R for 
periods up to 30 minutes. 
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The combustor and hydrogen-fuel injector (fig. 3) were both water- 
cooled. The injector was machined from a solid copper plate and had a 
blockage of 33 percent. Fuel was injected through 148 holes of 0.040- 
inch diameter facing downstream and drilled at 45° to the airflow direc- 
tion. Fuel was taken from a bank of high-pressure cylinders, reduced in 
pressure to about 150 pounds per square inch gage, and metered through a 
standard ASME orifice. Flow was controlled by throttling downstream of 
the metering orifice. 

The water-cooled nozzle had conical convergent and divergent sec- 
tions joined by a circular arc at the throat and an exit-to-throat area 
ratio of 9.5:1. The contour and instrument locations are shown in fig- 
ure 4. 


Instrumentation 

Static-pressure taps were drilled normal to the nozzle wall and were 
located spirally around the axis. Pressure taps were also located at the 
combustor inlet, outlet, and at the nozzle base. The pressures were read 
with the laboratory digital multiple-pressure recorder (ref. 17). Dif- 
ferential pressures at the metering orifices were read directly for air 
and through a mechanical pressure transmitter for the fuel. Water-orifice 
differential pressures and the fuel upstream-orifice pressure were meas- 
ured with transducers. After each change in fuel flow, the flow was al- 
lowed to reach steady- state conditions before any data were taken by 
monitoring the nozzle statics-pressures on a mercury manometer board. 

Temperatures of the cooling water, nozzle walls, air, and fuel at 
the metering orifices were measured with iron-constantan thermocouples. 

Air temperature at the combustor inlet was measured with a platinum - 
platinum-13 1 ^ rhodium cooled aspirated thermocouple probe. Two types were 
used during the course of this investigation: The first had a sonic noz- 

zle insert (ref. 18); and the second, which was a stagnation cup with 
bleed holes for the aspirated gas, proved to have a longer lifetime in 
the high-temperature gas stream since the junction was not exposed to a 
high-velocity stream. Radiation and conduction errors were negligible 
for both probes. A recovery correction of 1.5 percent was applied to the 
sonic aspirated probe. All thermocouple outputs were read with a digital 
voltmeter in conjunction with the laboratory central data system (ref. 

17). The inlet air temperature was also recorded on a strip-chart re- 
corder. 

Line reversal pyrometer . - The steady-state nozzle static tempera- 
tures at the port positions shown in figure 4 and the combustor-exit 
temperature profiles were measured with a self-balancing line reversal 
pyrometer (refs. 19 and 20) with a carbon-arc comparison source. Varia- 
tion of the source intensity was accomplished by rotation of an optical 
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wedge in the light path. The optical diagram for the instrument is 
shown in figure 5. The calibration procedure is given in appendix B. 

Both the D lines of sodium centered at 5893 A. and the blue cesium line 
at 4553 A were used for the reversal temperature measurements. The 
brightness temperature of the arc is higher at the lower wavelength and 
extended the range of the instrument to a higher temperature. Measure- 
ments could then be extended to the nozzle-iilet position. Sodium car- 
bonate or cesium sulfate was introduced into the test section as a powder 
in a carrier gas stream through a water-cool ;d probe. The probe was 
mounted in an actuator at the combustor exit and could be traversed radi- 
ally. The powder could thus be localized to a cross-sectional area of 
about 3/4-inch diameter at the nozzle inlet so that temperature profiles 
could be determined. The cesium sulfate used was somewhat hygroscopic 
(probably because of carbonate impurities) , so that care was required in 
preventing moisture from entering the powder feeder (ref. 20). The pow- 
der was dried before use, and 1 percent of Celite was added. Cylinder 
nitrogen was passed through a dessicant cartridge and used as the carrier 
gas for the powder. Air was used with the sodium carbonate. The nozzle 
ports were closed with contoured plugs except at the station where the 
optical unit was mounted. Lenses at the nozzle ports were used to reduce 
the field of view. 


Gas Sampling 


Combustion-gas samples were withdrawn at a distance 15 inches down- 
stream of the fuel injector and at various ladial locations through a 
0.020- inch-diameter water-cooled probe. A jressure ratio sufficient to 
obtain critical flow was maintained across the probe, and the samples 
were collected in 18-milliliter bottles. 

Analysis for hydrogen, nitrogen, and ojygen was performed with a 
Perkin-Elmer gas chromatograph using argon as the carrier gas. The con- 
centrations were determined by comparing pe£.k heights with known samples. 


Performance Calculat:' ons 

Equilibrium flow. - The thermodynamic state following constant- 
pressure combustion or isentropic expansion of air to an assigned pres- 
sure was determined by a set of equations defining the equilibrium, con- 
servation of mass, conservation of energy, partial pressures, and entropy. 
The successive approximation method was the same as that of reference 21. 
The thermodynamic data were taken from references 21 and 22. The calcu- 
lations were programed for the UNIVAC 1103 computer. 
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The temperature and composition following the combustion process are 
found for the measured initial conditions of air and fuel temperature , 
pressure, fuel-air ratio, and heat loss to the fuel-injector cooling 
water. The entropy following the isentropic expansion must be equal to 
the entropy before the expansion: 

S o = i Z["i( S T)i - Rn i In Pl J o 
i 

= [yi ( s § ) ± - Rn i ln p±J k d) 

i 

(All symbols are defined in appendix A.) If A is taken as the number 
of formula weights divided by the static pressure, 

Pi = n ± (2) 

Equation (l) can then be written in terms of partial pressures alone. 

The solution is found by successive adjustment of the temperatures and 
partial pressures to satisfy equation (l) for an assigned static pressure 
lower than the initial pressure. The velocity at any point is determined 
from the energy equation, 


u k = ^2gj(h 0 - h k ) (3) 

where h Q is taken at the nozzle inlet. The nozzle area at any step of 
the expansion is determined from the continuity equation 


_ / RT \ 

w UPbl 


(4) 


The throat area is determined by choosing the minimum from a pressure- 
area plot or simply as the minimum area determined from the stepwise cal- 
culation, using an initially small pressure step size. 

Frozen flow . - The composition is assumed frozen at the nozzle inlet 
following equilibrium constant-pressure combustion. The partial pressures 
at the k^ a step of the expansion are 

Pk 

(Pi)k = p^ (Pi^o 


( 5 ) 
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and the total entropy per equivalent formula weight of reactant is 

Pk / x 

s k = s 0 - R In — (6) 

Do 

where S G is given by equation (l). The temperature is determined by 

satisfying (l) for temperature, where the n-j ' s are the initial values 

following combustion. t?d 

Partially frozen flow . - The equilibrium calculations are followed 
to the step corresponding to the assumed freezing point. The composition 
is then assumed fixed, and the calculation is resumed as in the case of 
frozen flow but using the values of mole fraction at the freezing point 
as initial conditions. The initial enthalpy in equation (3) must con- 
tinue to be the value at the combustor. 


RESULTS AND DISCUSSICN 
Nozzle Calibration 

A cold-flow static-pressure calibration of the nozzle was made to 
determine the equivalent area ratio of the nczzle for the one -dimensional 
expansion of an ideal gas with y = 1.4. Dry air was used during the 
test. The boundary-layer displacement thickness was computed by the 
method of reference 23, described in appendi> C. Calculated displacement 
thicknesses for cold air, hot air, and stoichiometric combustion products 
are shown in figure 6. The equivalent area ratios from the cold-flow 
calibration and the area ratios corrected for the cold-air boundary -layer 
displacement thickness are shown in figure 7. The symbols show the area 
ratios at the static tap locations representing source flow and are the 
ratio of the spherical cap area to the sonic flow area. The sonic flow 
area was computed from the measured cold-air mass flow, total pressure, 
and temperature. Using this value for the throat area, the comparison 
to the calibration data corrected for boundary layer is quite good. 

Air-temperature measurements . - Radial air-temperature profiles were 
measured at the combustor outlet with the line reversal method using the 
tungsten lamp source and sodium tracer. The temperature was found to be 
flat within 120° F to 1 inch from the wall at about 3200° R. The average 
temperature was within 30° F of the thermocouple measurement. A fluctua- 
tion at a frequency of 3 to 4 cycles per second in the air- temperature 
measurements was sensed by the aspirating thermocouple. The magnitude of 
the fluctuation was of the order of ±100° F near the wall but was consid- 
erably reduced near the center of the duct. The magnitude of the fluc- 
tuations was also considerably reduced at the lower inlet pressure of 1.7 
atmospheres . 


9 ^ 21 - 
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Expansion of heated air . - A static-pressure profile for air at 
3300° R is shown in figure 8, together with the calculated equilibrium 
and frozen vibration profiles. The inlet temperature has been corrected 
for the cooling loss to the fuel injector. The equilibrium and frozen 
pressure ratios are evaluated at the equivalent one-dimensional area 
ratios. The area ratios at the pressure tap locations have been cor- 
rected for the change in boundary-layer displacement thickness from the 
cold- flow to heated-air conditions. The dotted line is plotted through 
the data by using the method of least squares and has a standard devia- 
tion about this line of 0.46 percent. It appears to indicate a lag in 
relaxation of vibrational energy. A discontinuity in the profile at an 
area ratio of 2.1 was typical of all the data including the calibration 
and was possibly due to a discontinuity in the nozzle contour. 


Combustion- gas samples . - A comparison of the measured and calcu- 
lated equilibrium composition of the samples plotted against sample 
equivalence ratio is shown in figure 9 and is presented in table I. 
Analysis was on a dry basis for hydrogen, nitrogen, and oxygen. Argon 
was used as the carrier gas and thus does not show up in the analysis, 
but was computed as 1.21 percent of the nitrogen present in the sample. 


Nitric oxide was assumed to be oxidized in the sample to NO^ and to 
be absorbed in the water present as the acids. It would presumably be 
lost by absorption in the desiccant used for drying the samples or as 
liquid remaining in the sample bottles. Corrections to the sample com- 
position for NO were made on the basis of the mole fraction present at 
equilibrium and the oxygen in the sample corrected for oxidation of NO 
to NOg* 


The water vapor formed during combustion was computed as 


n HgO = 2 An, 


Oc 


where the initial number of moles of oxygen is given as 


n r 


^2,0 3^6 p 2 , S 2 n N0 eq 

and the final number of moles as 


1 

+ v: niy 




n °2 = n 0 2 ,s + “ n ‘ 


1 

2 n N0, 


eq 


(7) 

( 8 ) 

(9) 


The equilibrium compositions were computed from the data of references 
21 and 22 at the same pressure, overall equivalence ratio, and combustion 
temperature as the sample. If the hydroxyl radical is assumed to recom- 
bine to form water, the compositions are in fairly good agreement and 
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indicate that there are no large deviations from equilibrium at the noz- 
zle inlet. 

The compositions are also given in table I at various radial sampling 
positions. It was not possible to hold conditions constant during a ra- 
dial survey because of the time involved in sampling and changing sample 
bottles; however, there are no pronounced discrepancies between the over- 
all and sample equivalence ratios at any radi^Q. position. 

td 

Combustor temperature profiles . - Radial temperature profiles at the 
combustor exit are shown in figure 10. These measurements were made by 
reversal of the 4553 A line of cesium while injecting cesium sulfate pow- 
der through a traversing probe. Superimposed on figure 10(a) is a sketch 
of the fuel injector showing its position relative to the powder-injecting 
probe. The asymmetry of the profile may be caused by a few plugged ori- 
fices and presumably would be averaged-out over other angular positions. 

Since the injector orifices were sized for critical flow at stoichiometric 
fuel flow at 3.6 atmospheres burner pressure, the distribution and mixing 
would be expected to be poorer at low fuel flews. This is evident from 
the profile shown for an equivalence ratio of 0.644 at 1.9 atmospheres 
in figure 10(b) . The theoretical combustion temperature considering only 
cooling loss to the fuel injector is noted on each profile and compares 
quite well with the measured temperature near the center of the duct for 
equivalence ratios near stoichiometric. The theoretical combustion tem- 
perature corrected for cooling loss to the combustor wall has been com- 
pared with the area- weighted average from the profiles and also agrees 
reasonably well for all the profiles. It appears, then, that the shape 
of the profiles for near- stoichiometric mixtures is due principally to 
heat transfer to the walls. 


Expansion of Hydrogen-Air Combusl ion Products 

Static -temperature measurement . - Sodium D-line reversal temperature 
measurements at an area ratio of 1.23 for pressures of 3.6 and 1.7 atmos- 
pheres and air temperatures from 2880° to 3188° R are plotted in figure 
11 as a function of equivalence ratio. The calculated equilibrium and 
frozen static temperatures are shown for comparison. The temperatures 
essentially follow the equilibrium curve at the higher pressure and show 
some departure from the equilibrium line at the lower pressure for close- 
to-stoichiometric mixtures. For lean mixtures the higher- than-equilibrium 
temperatures can be attributed to weighting by the reversal instrument. 
Since temperature zones higher than the calculated mean combustion tem- 
peratures exist because of poor fuel distribution and mixing, spreading 
of the injected tracer or sodium salts contaminating the air stream will 
tend to weight the indicated reversal strongly toward the high-temperature 
peaks. Another cause of this weighting is the effect of the air- 
temperature fluctuations noted earlier. Changes in mass flow to choke 
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the nozzle accompanying the air-temperature fluctuations produce changes 
in fuel-air ratio. This either accentuates or diminishes the combustion- 
temperature changes , depending on whether the mixture is leaner or richer 
than stoichiometric. The earlier results of this work reported in refer- 
ence 14 were more subject to error from this source, since the frequency 
response of the instrument used was limited by a 0.5-second lag of the 
tungsten filament. The frequency response of the modified reversal in- 
strument was limited to about 4 cycles per second by the servomotor that 
drives the optical wedge. This, however, was fast enough to respond to 
the changes just noted, so that averaging over the recorder fluctuations 
yielded the correct average temperature . Further errors at lean mixtures 
may be due to small-scale inhomogeneities in mixing that could not be re- 
solved because of the limiting time response or the spreading of the so- 
dium tracer stream. 

Figure 12 shows the measured reversal temperatures at an area ratio 
of 1.77. Here, the deviation from the equilibrium curve is quite pro- 
nounced for near-stoichiometric mixtures. The data at the lower pressure 
of 1.7 atmospheres indicate somewhat greater departure. The weighting 
toward high temperatures at lean equivalence ratios is also apparent at 
this location. 

The static temperatures are plotted in figure 13 as a function of 
area ratio at an equivalence ratio of 1. The area ratios for the data 
points were negligibly affected by the boundary- layer corrections and 
are taken from the cold-flow calibration at the axis of the optical path. 
Static temperatures at the last measuring station, where A/A* = 5.14, 
could not be measured because they were below the range of the instrument 
(2955° R) . The departure from the equilibrium curve occurs between the 
first and second measuring stations at 3.6 atmospheres, and somewhat be- 
fore the first measuring station at 1.7 atmospheres. 

Pressure measurements . - Static-pressure profiles for hydrogen-air 
combustion products are shown in figure 14(a) at an equivalence ratio of 
1.04, combustion temperature of 5230° R, and a total pressure of 3.6 at- 
mospheres. The area ratios at the static tap locations have been cor- 
rected for change in boundary-layer displacement thickness for a frozen 
boundary layer at an equivalence ratio of about 1.0. A curve has been 
fitted to the data by the method of least squares. The data have a 
standard deviation about this curve of 0.42 percent and indicate depar- 
ture from the equilibrium line slightly downstream of the throat. 

The pressure ratios at the last static tap location are shown plotted 
in figure 14(b) as a function of equivalence ratio. The frozen line has 
been corrected for boundary- layer losses and appears as the dashed line. 
The boundary layer was computed at equivalence ratios of 1.0, 0.6, and 
for heated air alone; and the correction was interpolated for intermediate 
values. Since the boundary- layer corrections depend very strongly on the 
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ratio of wall to free-stream enthalpies, the corrections would be sensi- 
tive to nonuniform free-stream properties. Ac lean equivalence ratios, 
where the profiles indicated higher temperatures nearer the walls due to 
uneven fuel distribution through the injector, the free-stream to wall 
temperature ratios higher than theoretical would result in a smaller dis- 
placement thickness; and hence the actual correction to the area ratio 
should be larger. The data would then follow the corrected frozen line 
fairly closely over the complete equivalence ratio range. 

M 

Calculation of the apparent freezing poin t. - Recent work at this 
laboratory (ref. 12) on gas sampling downstream of lean, flat, hydrogen- 
air flames has yielded a rate equation for an overall reaction: 

kp 

h 2 + | o 2 H 2 0 (10) 

k D 

A rate equation for the decay of hydrogen downstream of a flat flame was 
found to be best described by 

jfu i / 8100 

- I fL = 1.7X10 10 [H 2 ] 3 / 2 [0 2 e" RT (ll) 

dt 

This expression probably represents a lower limit, since diffusion was 
neglected in fitting the data. Since the concentrations of hydrogen 
downstream of these flames are in large exces ; of equilibrium, the rate 
equation is essentially a recombination rate and may be used to satisfy 
the Bray criterion (ref. 3) for a freezing point. The Bray freezing 
point is defined for the dissociation and three-body recombination of a 
diatomic molecule, where the net rate of reaction at equilibrium is ap- 
proximately equal to the rate of dissociation. Since r ^ = r R , where 
departures from equilibrium begin to occur, the Bray criterion for the 
freezing point may be written as 


- 2 * = Kr R (12) 

dt R 

where K is about unity. 

At equilibrium, the net rate of change of any species may be written 
as 

d[ ] ± d[ 1-j^eq d[ ] ± f( p) 
dt 


dt 


dt 


( 13 ) 
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where the concentrations are given by 


[ 



(14) 


Equation (13) states that the net rate due to reaction alone equals the 
total rate minus the change in partial molar density due to the expan- 
sion. The concentration due to density change will be 


1 1 ^p) - (15) 

M o 

The density is taken from the equilibrium flow calculations, since the 
reaction will affect the flow appreciably. A graphical solution of 
equations (ll) and (13) should give the freezing point. If the overall 
reaction given by equation (10) is considered, the net rate of change of 
[HgO] is determined from the equilibrium calculation as 

d[ h 2°] d f 1 

dt = U to { [H 2°]eq - [H 2 0] o f(p)j (l6) 

Equations (ll) and (16) are plotted in figure 15 against distance for 
stoichiometric mixtures at 3.6 and 1.7 atmospheres. A freezing point is 
indicated at an area ratio of 1.02 downstream of the throat at 3.6 atmos- 
pheres and at an area ratio of 1.07 upstream of the throat for a pressure 
of 1.7 atmospheres. At stagnation pressures of 3.6 and 1.7 atmospheres, 
freezing points were assumed at intervals In the vicinity of the throat, 
and the flow was calculated with frozen concentration assumed downstream 
of the point. The best fit of these calculations to the data was a 
freezing point at an area ratio of 1.1 for the higher pressure, and at 
an area ratio of 1.05 for the lower pressure. The computed temperatures 
are shown as dashed lines in figure 13. The computed pressures are 
shown by the square symbols in figure 14(a). Relating these apparent 
freezing points back to the rate equation (ll) yields a rate constant of 
2.3X1010 for the higher pressure and 6.4X10-^-^ for the lower pressure. 

In view of the uncertainty in determining the rate equation in the work 
on flat flames, and in extrapolating it to the nozzle flow conditions, 
the agreement is probably better than should be expected. The pressure 
data (fig. 14(a)) show good agreement with the calculations based on 
Bray’s approximation. For a more adequate comparison with the data, how- 
ever, an exact calculation of the flow with reaction is needed. 
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SUMMARY OF RESULTS 

Departures from chemical equilibrium have been observed for the ex- 
pansion of near-stoichiometric hydrogen-air conbustion products through 
a convergent-divergent exhaust nozzle. 

1. Static -temperature measurements show a nonequilibrium condition 
near the throat followed by a rapid decrease in temperature toward the 
frozen composition line. These measurements at 3.6 and 1.7 atmospheres 
indicated a greater degree of recombination at the higher total pressure. 

2. Wall static-pressure measurements also indicate a nonequilibrium 
condition slightly downstream of the throat for a total pressure of 3.6 
atmospheres; it was shown that the data approximated an expansion corre- 
sponding to one in which the composition had frozen at an area ratio of 
1.1 downstream of the throat. 

3. The application of an approximate freezing criterion to an over- 
all reaction rate for the oxidation of hydrogel and air predicts freezing 
of the composition in the vicinity of the nozzle throat. 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, June 13, 1961 
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APPENDIX A 
SYMBOLS 

A cross-sectional area of nozzle, number of equivalent formula 

weights 

a velocity of sound 

Ci, C 2 Planck T s radiation constants 
f function 

g gravitational constant 

H form factor, 6 */0 

incompressible form factor 
fp flat-plate incompressible form factor =1.3 
h enthalpy per unit mass 

J mechanical equivalent of heat 

K constant in eq. (12) 

k reaction-rate constant 

M Mach number 

jt mean molecular weight 

monochromatic areal intensity 
n number of moles 

Pr Prandtl number 

p static pressure 

p^ partial pressure of the i^* 1 species 

R universal gas constant 

r reaction rate 
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S entropy per mole 

T temperature 

T average temperature 

t time 

U velocity 

w mass -flow rate 

X mole fraction 

x axial distance from nozzle inlet 

y distance normal to axis of nozzle 

a fraction dissociated 

X ratio of specific heats 

6* displacement thickness 

€ emissivity 

9 momentum thickness 

Tv wavelength 

v kinematic viscosity 

p density 

t optical transmission factor 

cp equivalence ratio 

Subscripts: 

C cooled 

D dissociation 

e free stream 

eq equilibrium 


-1246 



E-1246 


15 


g combustion gas 

i species 

k k^* 1 step in expansion 

o stagnation state, initial conditions 

R recombination 

ref reference 

s sample 

T temperature 

tr transformed quantity 

w wall 

Superscripts: 

0 standard state 

* conditions at sonic throat 
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APPENDIX B 


CALIBRATION OF PYROMETER WITH ARC SOURCE 


By Donald Buchele 


The illumination system of the pyrometer with a carbon arc source 
is shown in figure 5. The tungsten lamp source of reference 19 has been 
replaced by a pair of circular absorption wedges, which are illuminated 
by a carbon arc image focused on them by a condensing lens. The logarithm 
of the wedge transmission factor varies from 0 to 3 and is linear with 
the angle of rotation. Wedge position is transmitted by a potentiometer. 
Self-balancing operation is obtained by rotation of one of the wedges. 

The other wedge is stationary and is used at the position of maximum 
transmission; it has an optical density increasing counterclockwise, 
opposite to that of the other wedge. This provides a large area of uni- 
form transmission for any position of the movable wedge. An interference 
filter and an air-cooled Plexiglas window are necessary to avoid exces- 
sive wedge temperature. Bandwidth of the interference filter is 3 per- 
cent of the wavelength. 




The intensity IK at the gas location and the corresponding tem- 


perature for the spectral wavelength used are determined by a calibration. 
Under blackbody conditions, the intensity of a radiating source is given 
by Planck's equation: 

N A,b = C R" 5 (e C ^ ATb - 1 (Bl) 


where T-^ is the temperature of the blackbody. For nonblackbody con- 
ditions, including an emissivity and an optical transmission factor 

between the radiating source and the poino of observation, the ob- 
served monochromatic areal intensity is 


N- 


A, a = T A e A N A,b 


= cy 


-5p/'T a . j)' 


(B2) 


where T a is an apparent temperature corresponding to the observed 


radiation intensity, 
cal ib rat ion. 


Equations (Bl) and (B2) are the basic equations of 


In the optical path of figure 5 radiation from the arc is trans- 


mitted by three groups of elements: 


X V r 2> an< ^ 


l 3' 


includes all 


optical elements except Tg, the interference filter; and is the 

nozzle or test-section window optics. Using nhe Wien approximation to 
Planck’s equation, N^ at the gas location is 

A -5 e ' C 2AT b = a _ 5 -C 2 /AT a 

\L A 2 A 3 1 


N A,a = C ^ 6 


l fc A T Ai T Ao T A. 


(B3) 
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where T b is the thermal temperature of the arc and T a is the apparent 
temperature of the arc observed at the gas location. An optical pyrom- 
eter sighted on the arc crater indicates a temperature T c at the pyrom- 
eter wavelength Ap and arc emissivity gp. The intensity is 


N. 


P; C 


= c l € p A p' 


-5 e - c 2/Vb = Cl XJ 5 e' C2/ Vc 


(B4) 


Letting = £p (the arc is assumed to be a gray-body radiator ), equa- 
tions (B3) and (B4) combine to yield 


1 1 . ( A P “ A ) , . A n 

T a T b C 2 ln € P “ C 2 ln T A, 1 T A,2 T A, 3 


(B5) 


This is the calibration equation. Wavelength A is of the spectral 
line; Ap is of the pyrometer, 0.65 micron. The emissivity of the arc 
is taken as 0.98 from reference 24. The transmission factors x^ p and 
T A, 3 vai T gradually with wavelength. They are measured photoelectrically 
in the vicinity of the spectral line by using only the light transmitted 
by the interference filter. The transmission factor x^ ^ of "the inter- 
ference filter is measured with a spectrometer using a bandwidth of 
0,0020 micron or less at the wavelength of the spectral line. A further 
correction in the case of the double-pass system for the transmission of 
the optics and reflector on the reflector side of the light path is made 
as in appendixes B and C of reference 18. The loop transmission factor 
resulting from reflection of light from the tungsten lamp filament would, 
however, be different in the case of the arc unit. A reflection factor 
of 0.12 is produced by the glass surfaces of the optical wedge, although 
the arc itself does not reflect. Consequently, a special value of xp 
used in the loop transmission term includes only the optics between the 
wedge and the gas. The correction for the loop transmission was found 
to be negligible for these measurements and was subsequently omitted. 
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APPENDIX C 


BOUNDARY- LAYER CORRECTIONS 

Momentum losses due to heat transfer to the cooled nozzle walls 
have Been calculated by using the method of reference 23, which was 
modified to include the effect of specific heat ratios other than 1.4. 
The method involves the momentum integral and moment -of -momentum equa- 
tions as simplified by using Stewartson's transformation. It assumes 
uniform and reversible free-stream properties with irreversible flow in 
the boundary layer. 


The quantity generally accepted as the test correction to apply to 
a flow is the boundary- layer displacement thickness 6*. It cannot be 
determined directly but may be calculated by using the relation 

H = ^ (Cl] 


where 6 is the boundary- layer momentum thickness and H is the form 
factor. These two quantities are determined with the transformed momen- 
tum and moment -of -momentum integral equations, respectively. Modified 
from reference 23, the equations are: 



t=d 
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and 



The quantities at the reference conditions are determined with the 
equation 



The enthalpies were computed with the data of reference 21. The calcula- 
tions were made with the assumption that the equilibrium composition at 
the combustor exit was frozen throughout the nozzle. Calculations for 
equilibrium composition throughout the nozzle were not made, because it 
was believed that the differences between a frozen boundary layer with a 
frozen free stream and one with an equilibrium free stream would be very 
small since Ng makes the largest contribution to the specific heat of 
the mixture. 

Physically, transformed quantities (subscript "tr") refer to a non- 
insulated boundary layer at a Mach number of zero and include the vari- 
ation of density with temperature. For boundary layers over insulated 
surfaces, and where Prandtl number equals 1, the transformed and incom- 
pressible (subscript ”i”) quantities are identical. 
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The compressible momentum thickness is obtained using the relation 


0 




Me 


Y~+l 

2 (r-i) 


'tr 



(C8) 


The relation between the compressible and incompressible form factors is 
as follows 


H = 



(% + l) + 




1 


(C9) 


A program was set up on the IBM 704 computer to solve equations 
(Cl) to (C9). The output from this program yielded the boundary-layer 
displacement thickness 5* as a function of x. The solution of equa- 
tion (C3) required the value of Hj_ at x = 0. This starting value of 
was determined by first evaluating the c cmpressible form factor H 

at Mach 0 using the tabulated integrals of reference 25. The incom- 
pressible form factor at Mach 0 was thei determined from equation 

(C9). Since the evaluation of in this manner included the effects 

of Prandtl number and heat transfer, it was relieved to be a better 
approximation than the use of the flat-plate value of 1.3. The values 
of Hj were of the order of 1.2 for both the heated air and combustion 
cases. 


The values of the transformed momentum thickness at x = 0 were 
calculated by the flat-plate analysis of reference 23. 

Referring to figure 6, note that a negative displacement thickness 
exists over the throat region for the heated air and combustion cases. 
Negative displacement thicknesses are to be expected where the thermal 
boundary layer is much thicker than the velocity layer, because higher 
densities in the boundary layer relative to those in the free stream re- 
sult in higher boundary- layer values of mass flow per unit area. Phys- 
ically, this fact means that the flow area is greater in this region 
than the calibrated area. 
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TABLE I. - COMBUSTION-GAS SAMPLES 


[p 0 , 3. 6 atm.J 



Combus- 

tion 

temper- 

ature, 

T, 

°R 

Equiv- 

alence 

ratio, 

9 

Radial 

probe 

posi- 

tion, 

in. 

Gas-composition mole 

fraction, 

% 2 

^2 

x h 2 

x h 2 o 

X ° H e, 

(eq)* 

5076 

1.16 


0. 605 

0.0034 

0.064 

0. 294 

0.0158 

(s)** 


1.09 

6 

. 618 

.0021 

.044 

. 326 


eq 

5197 

1.075 


. 614 

.0076 

.052 

.280 

.0224 

s 


1. 14 

6 

.611 

.0038 

.054 

. 318 


eq 

4694 

.577 


. 691 

.0684 

.0040 

. 200 

.0158 

s 


.666 

°2 

.674 

.0586 

.0004 

.243 


eq 

5237 

1.16 


.599 

.0050 

.0697 

.280 

.0211 

s 


.870 


. 660 

.0269 

.0009 

. 307 


eq 

5236 

1.02 

"1 

.616 

.0096 

.0491 

. 274 

.0245 

s 


.966 

% 

. 636 

.0164 

.0226 

. 308 


eq 

5168 

.915 


.637 

.0182 

.030 

. 263 

.0259 

s 


1.01 

0 

. 627 

.0110 

.030 

.312 


eq 

5188 

.945 


.632 

.0157 

.034 

. 266 

.0257 

s 


1.01 

2 

. 626 

.0157 

.034 

. 306 


eq 

5129 

.875 


. 644 

.0220 

.0246 

. 259 

.0255 

s 


.985 

1 

.631 

.0124 

.0273 

. 311 


eq 

5144 

.891 


. 641 

.0204 

. 0267 

.261 

.0257 

s 


.920 

2 

.641 

.021 

.0172 

. 301 



^Temperatures and composition calculated for equilibrium combustion 
at overall equivalence ratio. 

Composition and equivalence ratio calculated from sample. 
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Figure 1. - Ramjet net jet thrust for stoichiometric hydrogen-air combustion 
■products . 
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Figure 2. - Nozzle test facility. 
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Carbon arc 



Figure 5. - Optical diagram of line reversal pyrometer. 
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Figure 6. - Boundary- layer displacement thickness 
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Sample equivalence ratio, cp s 

Figure 9. * Gas sample composition, comparison wit a equilibrium conditions. Initial 
pressure, 3.6 atmospheres. 
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Reversal temperature 







Equivalence ratio 

(a) Initial pressure, 3.6 atmospheres; air temperature, 
2880° to 3185° R. 

Figure 11. - Reversal temperature measurements at area ratio 
of 1,23. 





Equivalence ratio 

(b) Initial pressure, 1.7 atmospheres; air temperature, 2885° to 2980° R. 

Figure 11. - Concluded. Reversal temperature measurements at area ratio 
of 1.23. 
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